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Abstract 

The sensitivity of THERA to different models of "new physics" has been stud- 
ied, both for the contact interaction approximation and for the resonance produc- 
tion. For contact interaction models conserving parity, scales up to about 18 TeV 
can be explored at THERA, extending considerably beyond the existing bounds. 
Significant improvement of existing limits is also expected for models with large 
extra dimensions. Effective Plank mass scales up to about 2.8 TeV can be probed. 
THERA will be the best machine to study leptoquark properties, for leptoquark 
masses up to about 1 TeV. It will be sensitive to the leptoquark Yukawa couplings 
down to Alq ~ 10~^. 

1 Introduction 

Search for "new physics" has always been one of the most important subjects in the field 
of particle physics. Possibility of discovering new particles, new interactions and/or other 
new phenomena is always considered as a main argument for building new, more powerful 
colliders. Collected in this paper are results concerning possible "new physics" searches at 
THERA, prepared as a part of the dedicated THERA physics study. The collider project 
and running options are briefly summarized in Sect. ^ Models considered in this paper 
are introduced in Sect. 0. The method used for the analysis has been developed for the 
global analysis of the existing data |1|, |^ and is briefly described in Sect. ^. Some of the 
results discussed in Sect. H have been already presented in M. 



2 THERA collider 



THERA has been proposed as the next generation ep collider, a straightforward extension 
of the TESLA project. It would bring high-energy electrons or positrons from the linear 



accelerator into collisions with high-energy protons from the existing HERA collider. 
Using both arms of TESLA center-of-mass energies of up to 1.6 TeV could be obtained 
at THERA, exceeding the energy range currently accessible at HERA by up to a factor 
of five. 

Following THERA running scenarios are considered in this paper: 

• For nominal electron beam energy of Ee=250 GeV and proton beam energy of 
Ep=1000 GeV integrated luminosity of about 40 pb~^ is expected in a year. Results 
presented for this option assume integrated luminosity of 100 pb~^ for e~p and/or 
100 pb~^ for e~^p collisions. It will be referred to as Thera-250. 

• Using both arms of TESLA, electron beam energy can be increased to Ee=500 GeV. 
Assumed integrated luminosity for this scenario is also 100 pb^^ for e^p and/or 100 
pb~^ for e~^p collisions. It will be referred to as Thera-500. 

• With TESLA machine upgraded in power, electron energies as high as i?e=800 GeV 
are possible for THERA operation. In this case proton beam energy is lowered to 
Ep=800 GeV, to provide maximum luminosity. Results presented for this option 
assume integrated luminosity of 200 pb^^ for e~p and/or 200 pb^^ for e'^p collisions. 
It will be referred to as Thera-800. 



3 Models of new physics 
3.1 Contact Interactions 

Four-fermion contact interactions are an effective theory, which allows us to describe, in 
the most general way, possible low energy effects coming from "new physics" at much 
higher energy scales. As very strong limits have been already placed on both scalar and 
tensor contact interactions 0], only vector contact interactions are considered in this 
analysis. The influence of the vector contact interactions on the ep NC DIS cross-section 
can be described as an additional term in the tree level eg — > eg scattering amplitude H, O : 
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where t = —Q^ is the Mandelstam variable describing the four-momentum transfer be- 
tween the electron and the quark, Cq is the electric charge of the quark in units of the 
elementary charge, the subscripts i and j label the chiralities of the initial lepton and 
quark respectively = L,R), and gj are electroweak couplings of the electron and 
the quark, and rj^j are the contact interaction couplings. 

In the most general case, vector contact interactions are described by 4 independent 
couplings for every quark flavor. As ep scattering is sensitive predominantly to electron- 
up and electron-down quark couplings, 8 independent couplings should be considered. 
However, it is not be possible, in one experiment, to put significant constraints on all 
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of these couplings simultaneously, without additional assumptions. Therefor, only one- 
parameter models, assuming fixed relations between the separate couplings, are considered 
in this paper. Relations between couplings assumed for different models are presented in 
Tab. |l| and Listed in Tab. |l| are models with defined coupling chirality. Models in 
Tab. ^ fulfill the relation 

eq I eq eq eq p. 

Vll + Vlr-Vrl-Vrr = 

which is imposed to conserve parity, and to avoid strong limits coming from atomic parity 
violation measurements. In the presented contact interaction analysis it is also assumed 
that all up type and down type quarks have the same contact interaction couplings: 

jjSd _ es _ eb 
'lij lij lij 

Coupling 7] can be related to the effective mass scale of contact interactions A: 
where the coupling strength of new interactions is by convention set to gci = v^ivr. 



3.2 Large Extra Dimensions 

Model proposed by Arkani-Hamed, Dimopoulos and Dvali 0, assumes the space-time is 
4 + n dimensional. Standard Model particles, including strong and electroweak bosons are 
confined to 4 dimensions, but the gravity can propagate in the extra dimensions as well. 
With very large extra dimensions, the effective Plank scale Ms can be of the order of TeV. 
The graviton, after summing the effects of its excitations in the extra dimensions, couples 
to the Standard Model particles with an effective strength of I/M5. At high energies, 
gravitation interaction can become comparable in strength to electroweak interactions. 
Virtual graviton exchange contribution to eq —>■ eq scattering can be described by an 
effective contact interactions. Contribution to the scattering amplitude (|l]), equivalent to 
the cross-section formula given in can be written as: 



eq 

Vll 


eq 

— Vrr — 


ttA , , , 


eq 

Vlr 


eq 

= Vrl = 





where t and s are the Mandelstam variables describing electron-quark scattering. By 
convention the coupling strength is set to A = ±1. 



3.3 Leptoquarks 

In this paper a general classification of leptoquark states proposed by Buchmiiller, Riickl 
and Wyler ||10| will be used. The Buchmiiller-Riickl-Wyler (BRW) model is based on the 
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Table 1: Relations between couplings for contact interaction models with defined coupling 
chirality considered in this paper. 
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Table 2: Relations between couplings for the parity conserving contact interaction models 
considered in this paper. 
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assumption that new interactions should respect the SU{3)c x SU{2)l x f/(l)y symmetry 
of the Standard Model. In addition leptoquark couplings are assumed to be family di- 
agonal (to avoid FCNC processes) and to conserve lepton and baryon numbers (to avoid 
rapid proton decay). Taking into account very strong bounds from rare decays [|ll|] it is 
also assumed that leptoquarks couple either to left- or to right-handed leptons. With all 
these assumptions there are 14 possible states (isospin singlets or multiplets) of scalar 
and vector leptoquarks. Table |^ lists these states according to the so-called Aachen no- 
tation An S(V) denotes a scalar (vector) leptoquark and the subscript denotes the 
weak isospin. When the leptoquark can couple to both right- and left-handed leptons, an 
additional superscript indicates the lepton chirality. A tilde is introduced to differenti- 
ate between leptoquarks with different hypercharge. Listed in Tab. ^ are the leptoquark 
fermion number F, electric charge Q, and the branching ratio to an electron-quark pair 
(or electron-antiquark pair), (3. The leptoquark branching fractions are predicted by the 
BRW model and are either 1, ^ or 0. For a given electron-quark branching ratio /?, the 
branching ratio to the neutrino-quark is by definition (1 — /5). Also included in Tab. ^ are 
the flavours and chiralities of the lepton-quark pairs coupling to a given leptoquark type. 
In three cases the squark flavours (in supersymmetric theories with broken R-parity) with 
corresponding couplings are also indicated. Present analysis takes into account only lepto- 
quarks which couple to the first-generation leptons (e, Ve) and first-generation quarks (m, 
d). It is also assumed that one of the leptoquark types gives the dominant contribution, 
as compared with other leptoquark states and that the interference between different lep- 
toquark states can be neglected. Using this simplifying assumption, different leptoquark 
types can be considered separately. Finally, it is assumed that different leptoquark states 
within isospin doublets and triplets have the same mass. 

In the limit of heavy leptoquark masses {M^q ^ -y/s) the effect of leptoquark produc- 
tion or exchange is equivalent to a vector type eeqq contact interaction. Contribution to 
the eg eq scattering amphtude (|I]) does not depend on the process kinematics and can 
be written as 



where Mlq is the leptoquark mass, Xlq the leptoquark-electron-quark Yukawa coupling 
and the coefficients a^J are given in Tab. | [|^. 

For leptoquark masses comparable with the available ep center-of-mass energy u- 
channel leptoquark exchange process and the s-channel leptoquark production have to 
be considered separately. Corresponding diagrams for F=0 and F=2 leptoquarks are 
shown in Fig. |l|. The leptoquark contribution to the scattering amplitude can be now 
described by the following formulae: 

• for M-channel leptoquark exchange ( F=0 leptoquark in e~q or e~^q scattering, or 
|F|=2 leptoquark in e^q or e~q scattering) 
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Table 3: A general classification of leptoquark states in the Buchmiiller-Riickl-Wyler 
model. Listed are the leptoquark fermion number, F, electric charge, Q (in units of 
elementary charge), the branching ratio to electron-quark (or electron- antiquark), (5 and 
the flavours of the coupled lepton-quark pairs. Also shown are possible squark assignments 
to the leptoquark states in the minimal supersymmetric theories with broken R-parity. 
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Table 4: Coefficients a^j defining the effective contact interaction couplings rj^J = a 
for different models of scalar (upper part of the table) and vector (lower part) leptoquarks. 
Empty places in the table correspond to alj = 0. 
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Figure 1: Diagrams describing leading order Standard Model processes and leptoquark 
contributions coming from F=0 and F=2 leptoquarks, for NC e^p DIS at THERA. 
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for s-channel leptoquark production {F=0 leptoquark in e'^q or e g scattering, or 
|-F|=2 leptoquark in e~q or e~^q scattering) 



where T is the total leptoquark width. The partial decay width for every decay 
channel is given by the formula: 

_ XIqMlq 

J- LQ 



87r(J + 2) ' 
where J is the leptoquark spin. 

For leptoquark masses smaller than the available ep center-of-mass energy direct pro- 
duction of single leptoquarks can be considered. In the narrow-width approximation, the 
cross-section for single F = 2 leptoquark production in electron-proton scattering (via 
the electron-quark fusion) is given by: 

a^P^^Q^{MLQ,XLQ) = {J + l)-^^-XLQq{xLQ,MlQ) (2) 
where q{x, Q"^) is the quark momentum distribution in the proton and Xlq = — 



4 Analysis method 

The analysis method used has been described in details in the recently published papers 
1^ . For all models considered, limits on the model parameters can be extracted from 
the measured distribution of NC DIS events at THERA. The leading-order doubly- 
differential cross-section for electron-proton NC DIS (e~p e~X) can be written as 
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where x is the Bjorken variable, describing the fraction of the proton momentum carried 
by the struck quark (antiquark), q{x,Q^) and qi^x^Q"^) are the quark and antiquark mo- 
mentum distribution functions in the proton and M^^ are the scattering amplitudes of 
([^), which can include contributions from "new physics". 

The cross-section integrated over the x and range of an experimental bin from 
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where i/max is an upper limit on the reconstructed Bjorken variable y, y = In the 
presented analysis this limit is set to ymax=0-95. Eq. ^ is used to calculate numbers of 
expected events in bins. Expected limits on model parameters, from high-Q^ NC e^p 
DIS at THERA, are calculated assuming that no deviations from the Standard Model 
predictions will be observed. 

For every value of the model parameter r] (±47r/A^ for the contact interaction models, 
±1/M| for large extra dimensions, {Xlq/Mlq)"^ for leptoquark models in the high-mass 
limitQ) the probability function describing the agreement between the model and the data 
is calculated: 

i 

The product runs over all Q"^ bins i (separately for e~p and e+p data). The probability 
Pi described by the Poisson distribution 

p.(^) _ ^(^)^ • exp{-n{r])) 

where and n{rj) are the measured and expected number of events in a given bin. This 
formula properly takes into account statistical errors in the measured event distributions. 
The systematic errors in the Standard Model expectations are assumed to be correlated 
to 100% between different Q"^ bins, and increase from 1% at Q^=1000 GeV^ to 5% at 
(5^=100000 GeV^. The method used to include systematic errors, as well as the migration 
corrections resulting from the assumed measurement resolution of 5% are discussed in 
detail in [0]. 

To constraint leptoquark Yukawa coupling values, for leptoquark masses smaller than 
the available ep center-of-mass energy, direct production of single leptoquarks is also 
considered, as described by (|^). Only the leptoquark signal in the electron-jet decay 
channel is taken into account. Expected signal from single leptoquark production, for 
given leptoquark mass Mlq and Yukawa coupling Xlq, is compared with the observed 
number of events from the Standard Model background (NC DIS) in the ±5% mass 
window. The background is suppressed by applying a cut on the Bjorken variable y, which 
is optimized for every leptoquark type as a function of the leptoquark mass. After the y cut 
is imposed, probability function V{Xlq, Miq) is calculated from the Poisson distribution 
(^. AsV is not a probability distribution, it does not satisfy any normalization condition. 
Instead it is convenient to rescale the probability function in such a way that for the 
Standard Model it has the value of 1: 

Vir] = 0) = 1. 

Using the probability function V{r]) limits on the model parameters are calculated. 
Rejected are all models (parameter values) which result in 

Viri) < 0.05 (5) 

^ For leptoquark masses comparable with the available center-of-mass energy, two parameter probabil- 
ity function V{Xlq, Mlq) is considered and limits on Xlq are calculated as a function of the leptoquark 
mass Mlq- 
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This is taken as the definition of the 95% confidence level (CL) exclusion limit.0 Exclusion 
limits presented in this paper are lower limits in case of mass scales A or M5, leptoquark 
mass Mlq or Mlq/Xlq, and upper limits in case of Xlq. For leptoquark masses smaller 
than the available center-of-mass energy, both indirect (from and direct Xlq limits 
are calculated, and the stronger one is presented. 



5 Results 



70 CL exclusion limits on the contact interaction mass scales A~ and (for negative 
and positive coupling signs) expected from the measurement of high-Q^ NC DIS cross- 
sections at THERA, are presented in Tab. |]and ^. Results presented are the mean values 
from 1000 MC experiments. Poisson fluctuations in the observed numbers of events can 
result in the statistical fluctuations in the limit values of the order of 10-20%. Current 
limits from the global contact interaction analysis fl^Q and from the global analysis of 



existing data in the large extra dimensions (ED) model |^ are included for comparison. 

For contact interaction models violating parity (models with deflned coupling chiral- 
ity; Tab. current limits from global analysis are already of the order of 10-20 TeV. 
This is mainly due to very strong constraints from the atomic parity violation (APV) 
measurements in cesium |16|, [1^. THERA running with the nominal electron beam 
energy of 250 GeV (THERA-250) will be only sensitive to mass scales from about 3 to 9 
TeV. With electron beam energy increased to 500 GeV (THERA-500), contact interaction 
mass scale limits improve on average by 20-30%. Another improvement by similar fac- 
tor is observed when going from THERA-500 to THERA-800 option. Nevertheless, even 
for high electron beam energies (THERA-500 and THERA-800) improvement of existing 
limits will only be possible for selected models (mainly models coupling to the u quark 
only). 

For contact interaction models conserving parity (Tab. |^), current limits from global 
analysis are, on average, lower than for parity violating models. At the same time THERA 
sensitivity increases. Already at THERA-250 mass scale limits can be improved for about 
half of the considered models, provided that both e~p and e~^p data are collected. With 
increasing electron beam energy, most limits can be signiflcantly improved. From com- 
bined e~p and e~^p data at THERA-800 limits on the contact interaction mass scales up 
to about 18 TeV can be obtained. 



^ For Gaussian shape of the probability function, condition corresponds to ±2.45fT limit. Mass 
scale limits presented in this paper would increase by 10 to 15%, if the definition more commonly used 
in the literature is used: V{r]) = 0.147 corresponding to ±1.96(t. However, this definition assumes the 
Gaussian shape of the probability function, which is not always the case. Therefor definition is used 
as more "conservative" . Same limit setting method has been used in the global analysis of existing data 

10- _ 

Numerical limit values presented in this paper differ slightly from limits presented in |]l4| . They have 

been recalculated using data on eeqq interactions only. Data from neutrino scattering experiments and 
from charged current processes, which can be included in the analysis when assuming SU{2)l x U{1)y 
symmetry of new interactions, were not used. This is because some of the considered models violate 
SU{2) invariance. 
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Table 5: 95% CL exclusion limits on the contact interaction mass scales A~ and A+ (for negative and positive coupling signs 
respectively) expected from the measurement of high-Q^ NC DIS cross-sections at THERA, for different running scenarios, as indicated 
in the table. Limits from the global analysis of existing data [1, 14] are included for comparison. 
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Table 6: 95% CL exclusion limits on the contact interaction mass scales A~ and A"*", and on the the effective Plank scales Mg and 
Mg in the large extra dimensions (ED) model (for negative and positive coupling signs respectively) expected from the measurement 
of high-Q^ NC DIS cross-sections at THERA, for different running scenarios, as indicated in the table. Limits from the global analysis 
of existing data [1, 14, 9] are included for comparison. 



For the model of large extra dimensions, THERA will improve existing limits in any 
configuration. This is because the graviton exchange contribution increases with the 
increasing center-of-mass energy. THERA-800 will be sensitive to the effective Plank 
scale Ms up to about 2.8 TeV. 

In the limit of heavy leptoquark masses {Mlq ^ ^/s) contact interaction model has 
been also used to set limits on the leptoquark mass to the coupling ratio Mlq/Xlq. 
Expected 95% CL exclusion limits, for different leptoquark models and different THERA 
running scenarios are presented in Tab. |^. Current limits from the global analysis [0 
are included for comparison. In most cases existing limits are already above THERA 
sensitivity, even in the highest electron energy option. Limits on Mlq/Xlq can be only 
improved for Vo and V1/2 models. 

For leptoquark masses comparable with the available center-of-mass energy, contact 



interaction approach has to be modified, as described in Sect. |3.3| . Limits on Xlq are 
calculated as a function of the leptoquark mass Mlq from the two-dimensional probability 
function V{Xlq, Mlq)- For leptoquark masses smaller than the available center-of-mass 
energy {Mlq < \/s), limits are also set from the measurement of the direct leptoquark 
production. It turns out both approaches give similar results [Q]. Measurement of the 
direct leptoquark production process results in better limits for low leptoquark masses 
{Mlq ^ a/s) and for leptoquark production involving valence quarks (production of F=2 
leptoquarks in e~p collisions or F=0 leptoquarks in e~^p collisions). measurement can 
results in slightly better limits (than expected from the direct production process) for 
leptoquark masses close to the center-of-mass energy and for leptoquark production from 
anti-quarks in the proton (F=0 leptoquarks in e~p or F=2 in e^p) . For leptoquark masses 
Mlq < \/s both kinds of limits are always calculated and the stronger one is taken. 

Shown in Fig. ||, ^ and ^ are expected 95% CL exclusion limits in {Xlq, Mlq), for 
different leptoquark models and THERA running with 250 GeV, 500 GeV and 800 GeV 
electron (positron) beam respectively. Limits expected from e~p data and from e^p data 
are compared. As expected, better limits on the F=2 leptoquark Yukawa coupling Xlq-, 
for Mlq < y/s, are obtained from e~p data, whereas e~^p data constrain better F=0 
leptoquarks. Differences between limits expected from e~p and e~^p data are smaller for 
scalar leptoquarks with Mlq > \/s. For high-mass vector leptoquarks it turns out that 
better limits can be obtained for "wrong" beam choice {e~p for F=0 leptoquarks and e'^p 
for F=2 leptoquarks). Limits expected from combined e~p and e^p data, for different 
THERA running scenarios, are compared with existing limits p| in Fig. ^ In all cases 
search for single leptoquark production at THERA significantly improves the existing 
limits. 

In Fig. 1^, 1^, H and ^, limits on the leptoquark Yukawa coupling Xlq and mass Mlq 
expected from THERA are compared with existing limits and limits expected from other 
future experiments 0, for 8^2^ Si, Vo and leptoquark models respectively.^ Lepto- 
quarks with masses up to about 2.0 TeV can be searched for at LHC, independently of 



Selected models were shown to describe the existing experimental data better than the Standard 
Model (|. 
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Table 7: 95% CL exclusion limits on M^q / X^q (in the limit of heavy leptoquark masses 
Mlq ^ V^) expected from the measurement of high-Q^ NC DIS cross-sections at 
THERA, for different running scenarios, as indicated in the table. 

Xlq- THERA will not be able to improve any limits if LHC excludes leptoquark masses 
below 1.6 TeV. However, if any leptoquark type state is discovered at LHC, THERA 
will be the best place to study its properties, covering the widest range in {Xlq, Mlq) 
space. Leptoquark mass, spin, fermion number and branching fraction (assuming lepto- 
quark decays into u + jet are reconstructed) can be determined. Yukawa coupling can be 
precisely measured down to the very small coupling values of the order of Xlq ~ 10~^, 
not accessible at LHC. 



6 Summary 

The sensitivity of THERA to different contact interaction models has been studied in 
detail. For models conserving parity, scales up to about 18 TeV can be probed at THERA, 
extending considerably beyond the existing bounds. Significant improvement of existing 
hmits is also expected for models with large extra dimensions. Effective Plank mass scales 
up to about 2.8 TeV can be probed. THERA will be the best machine to study leptoquark 
properties, for leptoquark masses up to about 1 TeV. It will be sensitive to the leptoquark 
Yukawa couplings down to Xlq ~ 10~^. 
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Figure 2: Expected 95% CL exclusion limits in {\lq, Mlq) space, for different leptoquark 
models (as indicated in the plot), for 100 pb~^ of e~p data (red curves) and 100 pb~^ of 
e+p data (blue curves) collected at 250 GeV electron (positron) beam energy (THERA- 
250). Limits based on single leptoquark production and high-Q^ NC DIS cross-section 
measurements. 
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Figure 3: Expected 95% CL exclusion limits in {Xlq, Mlq) space, for different leptoquark 
models (as indicated in the plot), for 100 pb~^ of e~p data (red curves) and 100 pb~^ of 
e~^p data (blue curves) collected at 500 GeV electron (positron) beam energy (THERA- 
500). Limits based on single leptoquark production and high-Q^ NC DIS cross-section 
measurements. 
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Figure 4: Expected 95% CL exclusion limits in {\lq, Mlq) space, for different leptoquark 
models (as indicated in the plot), for 200 pb~^ of e~p data (red curves) and 200 pb~^ of 
e~^p data (blue curves) collected at 800 GeV electron (positron) beam energy (THERA- 
800). Limits based on single leptoquark production and high-Q^ NC DIS cross-section 
measurements. 
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Figure 5: Expected 95% CL exclusion limits in {\lq, Mlq) space, for different leptoquark 
models and for different THERA running scenarios (as indicated in the plot). Limits 
based on single leptoquark production and high-Q^ NC DIS cross-section measurements 
from the combined e~p and e'^p data. Indicated in yellow are existing limits from global 
analysis 0]. 
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Figure 6: Comparison of expected 95% CL exclusion limits in {Xlq, Mlq) for 
leptoquark model, for different THERA running scenarios and other future experiments, 
as indicated in the plot. Presented limits correspond to 2x400 pb~^ of e^p data at HERA 
(yi=318 GeV), 2x100 pb'^ or 2x200 pb-^ of e^p data at THERA (7^=1.0,1.4 and 1.6 
TeV), 10 fb^^ of pp data at the Tevatron {^/s=2 TeV), 100 fb~^ of pp data at the LHC 
{y^=U TeV) and 100 fb'^ of e+e", 67 and 77 data at TESLA (^^=500 GeV). Also 
indicated are 95% CL exclusion limits from global analysis of existing data0. 
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Figure 7: Comparison of expected 95% CL exclusion limits in {X^q, M^q) for 5*1 lep- 
toquark model, for different THERA running scenarios and other future experiments, as 
indicated in the plot. Presented limits correspond to 2x400 pb~^ of e^p data at HERA 
(v^=318 GeV), 2x100 pb'^ or 2x200 pb-^ of e^p data at THERA (v/i=1.0,1.4 and 1.6 
TeV), 10 fb^^ of pp data at the Tevatron {^/s=2 TeV), 100 fb~^ of pp data at the LHC 
(^5=14 TeV) and 100 fb'^ of e+e", 67 and 77 data at TESLA (^^=500 GeV). Also 
indicated are 95% CL exclusion limits from global analysis of existing datap[. 
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Figure 8: Comparison of expected 95% CL exclusion limits in {X^q, M^q) for Vo lep- 
toquark model, for different THERA running scenarios and other future experiments, as 
indicated in the plot. Presented limits correspond to 2x400 pb~^ of e^p data at HERA 
(v^=318 GeV), 2x100 pb"^ or 2x200 pb-^ of e^p data at THERA (7^=1.0,1.4 and 1.6 
TeV), 10 fb^^ of pp data at the Tevatron {^/s=2 TeV), 100 fb~^ of pp data at the LHC 
{y^=U TeV) and 100 fb'^ of e+e", 67 and 77 data at TESLA (^^=500 GeV). Also 
indicated are 95% CL exclusion limits from global analysis of existing data0. 
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Mlq [GeV] 

Figure 9: Comparison of expected 95% CL exclusion limits in {X^q, Mlq) for 
leptoquark model, for different THERA running scenarios and other future experiments, 
as indicated in the plot. Presented limits correspond to 2x400 pb~^ of e^p data at HERA 
(yi=318 GeV), 2x100 pb'^ or 2x200 pb'^ of e^p data at THERA (v/i=1.0,1.4 and 1.6 
TeV), 10 fb^^ of pp data at the Tevatron {^/s=2 TeV), 100 fb~^ of pp data at the LHC 
(^5=14 TeV) and 100 fb'^ of e+e", 67 and 77 data at TESLA (^^=500 GeV). Also 
indicated are 95% CL exclusion limits from global analysis of existing datap[. 
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